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Abstract 

We show how to calculate correlation functions of two matrix models. Our method 
consists in making full use of the integrable hierarchies and their reductions, which 
were shown in previous papers to naturally appear in multi-matrix models. The 
second ingredient we use are the 14^-constraints. In fact an explicit solution of the 
relevant hierarchy, satisfying the M^-constraints (string equation), underlies the ex- 
plicit calculation of the correlation functions. In the course of our derivation we do 
not use any continuum limit tecnique. This allows us to find many solutions which 
are invisible to the latter technique. 
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1 Introduction 



Matrix models are believed to provide a (discrete) description of two dimensional 
gravity coupled to matter. Multi-matrix models might also provide a new approach 
to non-perturbative QCD. An accurate study of these models is therefore strongly 
motivated. While one-matrix models have been widely investigated, our knowledge 
of multi-matrix models is not as satisfactory (see however j^], 

11,111,111). 



In this paper we bring to a conclusion an analysis started in previous papers 
||],[||],0,[|l|] (see also ||], @). For simplicity we concentrate on two-matrix 
models and show how to find exact genus by genus solutions in these models Q By 
exact solutions we mean that we calculate the correlation functions of these theories 
without any approximation technique (except for genus expansion). The distinctive 
features of our approach are both its effectiveness - we are able to compute explicit 
solutions - and the method itself - we do not use any continuum limit tecnique. This 
allows us to find many solutions which are invisible to the latter technique, as will 
be clear in a moment. 

The essential tools used in this paper are the same introduced in the previous 
ones. In particular here we rely on the basic analysis carried out in |jl6|, although 
the emphasis on particular aspects may have changed in the meantime. 

For a model oi N x N hermitean matrices with bilinear coupling, we proceed 
as follows. 1) We associate to the initial matrix model partition function Zjy a 
discrete linear system (there are several equivalent ones), we identify the coupling 
conditions (string equations) and finally we specify how we can reconstruct Z^- from 
the quantities appearing in the linear system; we obtain in this way an extended 
description of our matrix model in the sense that this description is equivalent to the 
original path integral whenever the latter exists, but it is valid also for values of the 
coupling parameters for which the path integral is not defined. 2) We derive the Woo 
constraints. 

At this point we are faced with an integrable hierarchy of differential-difference 
equations; the solutions we are interested in are those that satisfy the above Woo 
constraints. To select them out we continue as follows. 

3) We separate the dependence on N from the dependence on the coupling pa- 



rameters via a procedure introduced in [16|, which consists in substituting the first 



flow equations into the remaining ones. This leaves us with a hierarchy of purely dif- 
ferential equations. 4) This new hierarchy is also integrable and contains integrable 
subsystems (reduced hierarchies) which can be nicely classified. 5) We identify the 
critical points. 6) We integrate the differential equations and plug the solutions into 
the Woo constraints to fix possible integration constants. In this way we obtain the 
sought for correlation functions. 

As one will see, all the above steps are of algorithmic nature. 

The results of our approach are of two types. The first, preliminary to the second, 
consists in the identification and analysis of the integrable hierarchy of differential 
equations which appear in two-matrix models. Since this is a long and self-consistent 



analysis we preferred to publish it in a separate paper [£l|. The second group of 



* Our approach is not of the type initiated by M.Kontsevitch although in a few cases the two 
approaches provide equivalent descriptions. 
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results is collected in this paper and concerns more specifically the calculus of cor- 
relation functions of two-matrix models. We are interested in the genus expansion 
of the solutions and concentrate in particular on genus results. They can be com- 
pared with the 2D gravity theory in the continuum and with topological field theory 
models. 

There are two broad families of models and correlation functions (CF's). To 
understand this point one should remember that two-matrix models describe the 
interaction of the first matrix Mi with itself, which is described by a (possibly infinite) 
polynomial in Mi, of the second matrix M2 with itself, which is described by a similar 
polynomial in M2, and finally by a bilinear interaction M1M2 with a coupling c. The 
first family of models is the one in which c 7^ 0, i.e. we consider the whole theory. 
The second family consists of models in which either c = or c does not appear, 
and the M2 sector is disregarded. We provide one example of the first family, but 
we mostly concentrate on the second. In the latter the matrix M2 plays a spectator 
role; its self-interaction and interaction with Mi are however essential. This goes as 
follows. Every model is characterized by an integrahle hierarchy, which is a reduction 
of the KP hierarchy. The order of the second potential determine a subclass of 
these models: if the order is p then the model is characterized by either a 2p-field 
representation of the KP hierarchy or by one of its integrable reductions. The models 
are accordingly labeled A^j, with I = 0, ... ,p — 1, where M.'^~^ corresponds to the 
just mentioned 2p-field representation of KP, while the others are reductions. In 
particular corresponds to the p-th KdV hierarchy. The latter are characterized 
by the fact that the correlation functions do not depend on the size of the matrices, 
A^. For all the other models the correlation functions do depend on A^, therefore 
these models cannot be 'seen' with continuum limit techniques. The first potential 
has another role, the non-vanishing coupling parameters in it determine the small 
phase space and the critical points. 

We remark that an explicit solution of the relevant hierarchy, satisfying the W- 
constraints (string equation), underlies the explicit calculation of the correlation 
functions of a given model. The solutions of the hierarchies relevant for CF's are 
characterized by a genus by genus homogeneity, defined according to a degree which 
is determined by the particular critical point one is considering. This is a distinctive 
feature of enormous practical utility. 

In elaborating our method we have been inspired by earlier works, in particular 
[^],[^],[^] and [^]. The models we obtain are very close to the topological field 
theory models coupled to gravity studied therein. 

This paper is organized as follows. In section 2 we review the relevant results of 
In section 3 we show how one can compute CF's directly from the VF-constraints; 
we show an example of CF's with c ^ 0. In section 4 we collect the relevant results of 
and [21|. Section 5 is entirely devoted to the explicit calculation of genus CF's 
at the first critical point via integration of the relevant hierarchy in some specific 
models, precisely the three models specified by the 4-field representation of the KP- 
hierarchy and its integrable reductions. In section 6 we give a general prescription 
for calculating CF's and a few more explicit examples. Section 7 is devoted to 
higher critical points: we show again how to calculate the corresponding correlation 
functions. The last section is devoted to a brief discussion of the connection of the 
previous models with topological field theories and 2D gravity coupled to matter. 
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The Appendices contain explicit expressions for VF-constraints, flow equations and 
CF's for the various models. 



2 Review of previous results 

For simplicity we limit ourselves to two-matrix models. They are initially defined by 
the partition function 

Zjv(t,c) = J dMidM2e^~'^ 
where Mi and M2 are Hermitian N x N matrices and 

U = Vi + V2 + CM1M2 

with potentials 

Pa 



K. = EWK a = 1,2. (2.1) 

r=l 



The paS are finite positive integers. 

A clarification is in order concerning the coupling constants ta^r- Later on we 
will let Pa ^ 00, thereby introducing an infinite number of couplings. Therefore the 



coupling constants split into two sets, those appearing in eq.(2.1) (which define the 
model), and the remaining ones (which are introduced only for computational pur- 
poses). In terms of ordinary field theory the former are analogous to the interaction 
couplings, while the latter correspond to external currents (coupled to composite 
operators). In the case of matrix models we will not make any formal distinction 
between them. This will allow us to obtain very symmetrical and powerful formulas 
- the W^-constraints for example. Of course the distinction is substantial and will 
appear when calculating the correlation functions. These will depend only on the 
interaction couplings, while the external ones will be set to zero once they have done 
their job. 

Let us return to the models defined at the beginning. Our purpose is to study 
the their correlation functions (CF's), i.e. the correlation functions of the operators 

TrMf, TrM|, Tr(MiM2) 

See the beginning of section 3 for more precise definitions. 

The ordinary procedure to calculate the partition function consists of three steps 
[^],[^],[^]: (i) one integrates out the angular parts so that only the integrations 
over the eigenvalues are left; {ii) one introduces the orthogonal polynomials 

^n.(Ai) = A" + lower powers, Vni^2) = A2 + lower powers 

which satisfy the orthogonality relations 

dAidA2en(Ai)e^i(^i)+^2(^2)+^^i^2r?„(A2) = K{t, c)5nm (2.2) 
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(iii), using the orthogonahty relation ( |2.2| ) and the properties of the Vandermonde 
determinants, one can easily calculate the partition function 

N-l 

ZN{t,c) = const Nl Yl hi (2.3) 

i=0 

Knowing the partition function means knowing the coefficients hn{t,cys. 

The information concerning the latter can be encoded in a suitable linear system 
plus some coupling conditions, together with a relation that allows us to identify Zjy. 
But before we pass to that we need some convenient notations. For any matrix M, 
we define 

iM)i~ M,j^, Mij = Mji, Mi{j) = M^j_i. 

As usual we introduce the natural gradation 

deg[Eij] =j-i 

and, for any given matrix M, if all its non-zero elements have degrees in the interval 
[a,b], then we will simply write: M € [a,b]. Moreover will denote the upper 
triangular part of M (including the main diagonal), while Af_ = M — M+. We will 
write 

N-l 

Tr(M) = J2 Ma 

Analogously, if L is a pseudodifferential operator, L+ means the purely differential 
part of it, while L_ = L — L+. 

Let us come now to the step 1) mentioned in the introduction. First it is con- 
venient to pass from the basis of orthogonal polynomials to the basis of orthogonal 
functions 

*n(Ai) = e^i(^^)c„(Ai), <^n{X2) = e^^(^^)77n(A2). 
The orthogonality relation becomes 

J dAidA2^'„(Ai)e^^i^2«>„(A2) = 5nnMt,c). (2.4) 

As usual we will denote the semi-infinite column vectors with components ^o, ^i; ^2, • • • , 
and ^o,^i,^2t ■ ■ , by ^' and respectively. 

Next we introduce the following Q~type matrices 

J dXidX2^n{Xl)Xae^^'^^^m{X2) = Qnm{a)hm = Qmn{a)K, a = 1,2. (2.5) 

Both Q{1) and Q(2) are Jacobi matrices: their pure upper triangular part is /+ = 
J2i Ei^i+i- 

Beside the above Q matrices, we will need two P~type matrices, defined by 

J dAldA2(^^„(Al))e^^^^^«>^(A2) = Pnm{l)hm (2.6) 
J dXidX2^n{Xl)e'^'^'( — <^>m{X2))=Pmn{2)hn (2.7) 
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The two matrices ( p.SD we introduced above are not completely independent. 
More precisely both (5(a) 's can be expressed in terms of only one of them and one 
matrix P. Expressing the trivial fact that the integral of the total derivative of the 
integrand in eq.(2.4) with respect to Ai and A2 vanishes, we can easily derive the 
constraints or coupling conditions 



P{1) + ci2Qi2) = 0, 



ci2Q(l) + ^(2) = 0, 



It is just these coupling conditions that lead to the famous Tyi+co~constraints on the 
partition function at the discrete level. We will also refer to them at times as string 
equations. From them it follows at once that 



Q{a) G [ 



where 



and 



mi = p2 



a 



7712 



1,2 



1, 



"2 = Pi - 1 



These equations show that a finite band structure for the Q{a) is only allowed if the 
number of perturbations is finite. Conversely, if we want, say, Q{1) to possess the 
full discrete KP structure - we should let p2 — > 00. 

The derivation of the linear systems associated to our matrix model is very simple. 
We take the derivatives of eqs.( |2.4| ) with respect to the time parameters ta,r, and use 
eqs.(p.5|). We get in this way the time evolution of ^' and ^, which can be represented 
in two different ways: 

(*) Discrete Linear System I: 



r Q(l)^'(Ai) = Ai'f (Ai), 



a 



dt 

a 



-Q^(2)M/(Ai 



l<k<pi, 
l<k<p2, 



[ |^vI/(Ai) = P(l)^'(Ai). 
The corresponding consistency conditions are 



^Q(1) = [Q(1), Qlia)] 
^ P(1) = [P(1), Qlia)] 



dta,. 



(2.9a) 
(2.9b) 

(2.9c) 



( 2.9b| ) and (2.9c) are discrete KP hierarchies, whose integrability and meaning were 



discussed in [17| 
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(**) Discrete Linear System II : 



Q(2)$(A2) = A2$(A2), 

at.^iM) = Q^(2)$(A2), 

^$(Ai) = -Q'i(l)$(A2), l<k<pi 



(2.10) 



_ g|^$(A2) = P(2)$(A2). 
with consistency conditions 



[Q(2), P(2)] = l, (2.11a) 
^ Q(2) = [Q(2), Q^(a)] (2.11b) 



d 

dta,k 



P(2) = [P(2), Qlia)] (2.11c) 



The third element announced in the introduction is the hnk between the quantities 
that appear in the hnear system and in the couphng conditions with the original 
partition function. We have 

lnZjv(i, c) = Tv(Q'-(a)), a = 1, 2 (2.12) 



da, 

It is evident that, by using the flow equations above we can express all the derivatives 
of Zn in terms of the elements of the Q matrices. For example 



lnZN{t,c) = (Q'{a))^^_^, (2.13) 



^<r<Pa; a = 1,2 
We also recall the coupling dependence of the partition function 

— lnZAr(t, c) = Tv(q(1)Q(2)) (2.14) 

Knowing all the derivatives with respect to the coupling parameters we can recon- 
struct the partition function up to an overall integration constant. 

We also remark that, since the RHS of the above equations is always defined, 
they give us a definition of Zjy even in subsets of the parameter space where the 
path-integral is ill-defined. 

We will be using the following coordinatization of the Jacobi matrices 

mi m2 

Q(i) = /+ + E E 2(2) = ^+ + E E bimi,i-i (2.15) 

i /=0 i 1=0 

One can immediately see that 
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As a consequence of this coordinatization, eq. (|2.13|) gives in particular the important 
relation 

lnZN{t,c) = ai{N), (2.16) 



To end this subsection we write down explicitly the and t2, inflows, which 
will play a very important role in what follows 

Q^^^iU) = + 1) - o/+i(j) + aKi)(ao(i) - ao{j - O) (2.17a) 

^ai{j) = Rj-i+iai^iU) - Rjai^iii - 1) (2.17b) 

j^biij) = Rj-i+ihi-i{j) - Rjbi^iij - 1) (2.17c) 

-^Hj) = h+iU + 1) - 6z+i(j) + bi{j){bo{j) - boU - I)) (2.17d) 



2.1 constraints 

The constraints (or simply VF-constraints) on the partition function for our 

two-matrix model were obtained in |l^ by putting together both coupling conditions 
and consistency conditions (see above). In other words the Wi+00 constraints contain 
all the available information. They take the form 

wMz^(t,c) = 0, #Mz^(t,c) = r>0; n > -r, (2.18) 

where 

T^M ^ (-c)"/:H(l)-£[:r^(2) (2.19a) 

H^M ^ {-crd:\2)-C^T\l) (2.19b) 



The generators Cn\l) are differential operators involving N and ti^k, while Cn\2) 



have the same form with ti^k replaced by t2,k- One of the remarkable aspects of ( ^.18 ) 
is that the dependence on the coupling c is nicely factorized. The Cn\l) satisfy the 
following 1^1+00 algebra 

[/:W(1),£W(1)] = (n - m)/:™ „(1) (2.20a) 
[/:[?](1),£W(1)] = (n - 2m)/:lfl.„(l) + m{m + 1)/:L'U(1) (2.20b) 
[/:[?](1),£P1(1)] = 2(n - m)/:L'U(l) - (n - m)(n + m + 3)/:L'i^(lX2.20c) 

and in general 

[d;;\l),d^{l)]={sn-rm)cttL~'\^) + ---, (2.21) 
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for r, s > 1; n> —r,m> —s. Here dots denote lower than r + s — 1 rank operators. 
We notice that this Vl^i+oo algebra is not simple, as it contains a Virasoro subalgebra 
spanned by the £n'(l)'s. We sec that once wc know these generators and £[^2(1)) 
the remaining ones are produced by the algebra itself. The explicit expression of the 
basic generators is given in Appendix Al. 

The algebra of the (2) is just a copy of the above one, and the algebra satisfied 
by the Wn and by WA is isomorphic to both. We refer to this abstract algebra with 
the symbol W. 



3 How to compute correlation functions from 
the VK— constraints 



This section is a sort of intermezzo. It does not represent our final solution of the 
problem. But still it will teach us a lot. We said above that the PF-constraints contain 
all the available information of our two-matrix model. Therefore, starting from them 
only, we should be able to calculate all the correlation functions of the model. This 
is indeed true and we are going to give in the following a sample calculation of this 
fact. To simplify the notation we set 



From the definition of the partition function we inherit the definition of the correlation 
functions 



These formulae do not make any distinction between internal and external couplings 
(see the specification at the beginning of section 2) . Eventually the external couplings 
will be set to zero. 

Now wc write the M^-constraints in this new language. wi^|Zjv = and H^l^jZjv = 
become, respectively 




TVM| = ak 



< Tfei • ■ ■ rfe„ » 




00 



^ ktk < Tk-i » +Nti + c < ri »= 



(3.1a) 



k=2 



00 



^ ksk < (Tk-i > +Nsi + c < (71 >= 

k=2 



(3.1b) 



Instead Wq 'Zjsf = and Wq 'Zn = give rise to the same equation 



00 



00 




(3.2) 
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The constraint M/'l^'z^v = becomes 



c(^/t«<r/+i»+(iV + l)<Ti»)+ J2 hhsi.si^ <^ ai,+i^_i :^ 
1=1 li,h=l 
oo 1-2 

+ (^km-k-i » + < (Tfc ai-k-1 » ) 

i=3 k=l 

oo 

+{2N + l)Y^ Isi < ai_i > +(iV2 + ^ (3.3a) 

1=2 

It should be clear by now that W^^Zn = gives an exactly symmetrical equation 
where tk and are interchanged with Sk and Ufc, respectively. In a similar way we 
can write all the other VF-constraints. 

Let us consider now the problem at issue in a simplified situation, i.e. in genus 0. 
This will be enough to shed light on the main features. In this case the W^-constraints 
come out simplified, p^. To see this assign a suitable degree to the couplings 



deg=[], [t^] = [sk]=x-k, [N]=x, [c]=x-2 (3.4) 
and the degree 

= 2x (3.5) 

to the genus zero part of the "free energy" (let us set from now on InZ = T). Here 
X is, for the time being, an unspecified positive real number. These assignments 
follow directly from the -constraints and the first flow equations. We have only 
arbitrarily (but without loss of generality) fixed the scale. We will be looking for 
CF's homogeneous (genus by genus) in the couplings in accordance with the above 
degree. 

In the ly-constraints the genus contribution is the homogeneous part of highest 
degree. In any correlation function we should henceforth append a label <C to 
indicate this contribution. However we will avoid this cumbersome operation by 
understanding that the correlation functions in the remaining part of this section 
refer only to the genus contribution. With this convention the genus version of 
( p.la|) , ( |3.1b| ) and ( |3.2D remain (formally) the same, while (3.3a) takes the form 



1=1 li,h=l 
oo 1—2 oo 

+ ^/sz^ < CTfe > +2iV^/sz < » +iV2si = (3.6) 

/=3 k=l 1=2 

and so on. 

The next step consists in specifying the critical point (or the model) we want 
to consider. This subject will be explained in detail later. For the time being we 
consider the first critical point of the simplest possible model: 

2t2 = -1 tfc = 0, k>2, 2s2 = -1 Sk = 0, k>2, (3.7) 
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Therefore CF's will be functions of N, c,ti = t and si = s only. To distinguish this 
from the general case we will denote the CF's with the symbol < • > instead of 
<^ ■ In order to preserve homogeneity, we set 



X 



in the above degree assignments. The genus -constraints ( 3.1a ), ( p.lbD 
and ( |3.6|) become respectively 

- < n > +Nt + c < cji >= 0, 

- < CJl > +Ns + c < Ti >= 0, 

- < T2 > +t < n >= - < a2 > +s < ai > , 

c(- <T3> +t <T2> +N < Ti >) 

+ < 0-3 > -2s < 0-2 > +(s^ - 2N) < cTi > +iV^s = 

Writing down the appropriate formulas for the other constraints one quickly realizes 
that these formulas form a recursive and overdetermined system of algebraic equa- 
tions for the one-point correlation functions. A simple computer program allow us 
to explicitly calculate them. For example we obtain 



<Ti>=N 



t + cs 



<ai>= N 



s + ct 



(3.9) 



and so on (the first few are written down in Appendix CI). We immediately remark 
that, setting c = and s = in these formulas, we obtain the CF's of the NLS model 



already met in the context of one-matrix model, \1G\- The NLS model correspond to 
A^2 ™ the classification of the present paper. 

One can also derive the multi-point correlation functions by simply differentiating 
an appropriate number of times the M^-constraints with respect to the couplings 
and repeating the same procedure. However this method of calculating correlation 
functions is not very economical since it relies on our being able to solve systems 
of algebraic equations which become more and more complicated the higher the 
critical point is. There is also another reason why this method is insufficient: since 
we do not have at our disposal a theory of the reduction of the W algebra and of 
the VF-constraints, it is impossible to study the reductions of our system, i.e. to 
study consistent subsystems which, as we shall see, constitute a large part of the rich 
structure of two matrix models. 

We need another method and this is what we are going to study in the following 
sections. However from the above example we have extracted some useful informa- 
tion: 

1) integrability of the discrete hierarchy on which the ly-constraints are based 
translates itself into a recursive set of equations for the CF's; 

2) setting c = (at the end of the computation) and disregarding and a^, 
reduces the CF's to those of a simpler model, the NLS model studied in jl^ - as will 
be seen, this is a general fact. 



3) it follows from |16| that, in the same conditions as at the previous point, the 
constraints (2.18) can be replaced by the much simpler ones 



LnZ 



(3.10) 
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n-1 q2 

+ Y.7^^ + N^Ko + Nti6n,-i, n>-l 

The -Ln's satisfies the Virasoro algebra and these Virasoro constraints are enough to 
completely determine the correlation functions. Effective W constraints of this type 
will appear also in other cases (other critical points, other models). 



4 Differential Hierarchies of Two— Matrix 
Models. 

The method alluded to at the end of the previous section is based on extracting 
differential hierarchies of flow equations from the discrete ones that characterize two- 
matrix models. Let us resume the summary of section 2. We recalled there that 
two-matrix models can be represented by means of coupled discrete linear systems, 
whose consistency conditions give rise to discrete KP hierarchies and string equations. 
Here we review the method introduced in [^] and used in [^] and [^], to transform 
the discrete linear systems into equivalent differential systems whose consistency 
conditions are purely differential hierarchies. This is tantamount to separating the 
dependence from the dependence on the couplings. 

The clue to the construction are the first flows, i.e. the ti^i and t2,i flows. For the 
sake of simplicity let us consider the system / and the flow ( p. 17a ). Let us consider 



the generic situation in which Q(l) has mi = p2 — ^ lower diagonal lines (see the 
parametrization ( 2.15| )). To begin with let us notice that 



dti 

and let us adopt for any function f[t) the convention 
We can rewrite 



d 

= ^n+i + ao{n)^n (4.1) 



= (4.2) 



where 



-1 CXJ 

^ 7^ = E(«o(n)5-^)' (4.3) 

Then it is an easy exercise to prove that the discrete spectral equation 

Q(l)^(Ai) = AiM/(Ai) 
is transformed into the pseudo differential one 

Ln{l)^n = Al^-n (4.4) 
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where 



mi 



1 1 



(4.5) 



1=1 

mi 



5 + Vai(n)- 

o — ao{n — I) o — aQ{n — I + 1) 



d — a(j{n — 1) 



Ln{l) is an operator of KP type |^. Actually it is in general a reduction of the KP 
operator 

oo 

LKP = d + J2wid-' 
1=1 

where wi are generic (i.e. unrestricted) coordinates. 

Proceeding in the same way for the other equations of system I we obtain the new 
system in differential form 



L„(l)^'„ = Ai^-^ 



d 



9*1, r 



(4.6) 



aAi 



In particular we have the following consistency conditions 

^ L„(1) = [(L;(1))+,L„(1)], 

L41) = [L„(1),(L;(2))_], 



d 



dt 



2,r 



(4.7) 
(4.8) 



They are purely differential equations. 

Let us come now to the n dependence of the above equations. The operator L„(l) 
introduced above (l4.5|) depends on the coordinates of many sectors. We showed in 
[0] that, using the first flow equations, we can transform it into a pseudo differential 
operator depending only on the n-th sector coordinates. This was implemented at 
the price of introducing a very complicated expression. For our purposes in this paper 
this is not very convenient and we had better proceed in another way. Precisely we 
introduce mi "fields" S*!, Sm^, related to the "field" oq in the following way 



Si{n) = ao{n - i) 
Then we can rewrite Ln(l) in the following way 

1 1 

1=1 



(4.9) 



mi 



Lnil) = d + Y,ai{n] 



d-Si{n) d-Si_i{n) d-Si{n) 



*In p7| we imposed the conditions a/(n) = 0, n < I, but actually there is no reason to impose these 
conditions on the basis of the matrix model, therefore we drop them here. Consequently the discussion 
about non-universality in section 6 of must be dropped too; in the language of that section, the only 
alternative present in multi-matrix model is the universal one. 
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So we have achieved the same result as m [|^, except that, of course, the fields Si 
are not independent. But following a long tradition in field theory, we will consider 
these fields as completely independent from one another in all the intermediate steps 
of our calculations and only eventually impose the condition ( [4.9D . 

To further simplify the notation we will consider henceforth the lattice label n on 
the same footing as the couplings and write 

ai{n, ...) = aj(n)(...), Si{n, ...) = Si{n){...) 

where dots denote the dependence on ti,fc,t2,fe and c. So the expression of Li{n) get 
further simplified to 

m -, , -, 

where, for simplicity, we have dropped the label (1) too. A similar simplification has 
to be understood also for the other equations of the system I above. This simplified 
form is the one we constantly refer to throughout the remaining part of the paper. 
In particular we recall the integrable hierarchy 



_d_ 

dtr 



L=[(r),,L], (4.11) 



and = ti^r- 

We can of course extract an analogous pseudo-differential operator for system II, 
by using the first flows ( |2.17d| ). One can also reconstruct the discrete systems from 
the differential ones and the first flows, for an explicit example see [pO[]. 



This hierarchy and the relevant integrable reductions have been studied in |18] 



and in the companion of this paper |21|. There we showed that to each operator 



such as ( [4.101 ) there correspond m distinct integrable reductions which are obtained 
via Dirac reduction procedure by successively suppressing the S fields one by one 
(the order being irrelevant). Of each of them we gave a Lax operator representation, 
which not only constitutes a very quick proof of integrability, but also provides a 
very efficient way to explicitly calculate all flows. To each such Lax operator, as we 
shall see, we can associate a full set of correlation functions, i.e. a model. Therefore 
we can say that for any p = p2 = m + 1 {p2 is the order of the potential V2) there 
correspond p distinct integrable models A^p, where I runs from p to and counts 
the number of S fields. The original unreduced model M.^~^ is defined of course by 
the Lax operator ( |4.10D above. The most reduced model is characterized by the 
p-th KdV Lax operator and the p-th KdV hierarchy. 

To end this section let us comment a bit on the naturalness of the reduction 
procedure. Our reduction procedure consists in restricting the second Hamiltonian 
according to the chosen constraints and calculating the second Poisson brackets ac- 
cording to the Dirac recipe. The rest is automatic. In the framework of the Hamil- 
tonian systems this reduction procedure is completely natural and allows us to select 
consistent subsystems of a given integrable system, or, in other words, to find solu- 
tions of the initial systems in which a part of the degrees of freedom are disregarded. 
We stressed in section 2 that our unreduced systems are equivalent to the initial path 
integral and that they even provide an enlarged definition of the latter. However we 
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do not know in general how to go back from a reduced system to the path integral 
formulation. We know that in the simplest case (one-matrix model giving rise to 
the NLS hierarchy, see [|l6|) the reduction is the KdV hierarchy and that the KdV 
hierarchy corresponds to the even potential case via double scaling limit. Perhaps a 
generalization of this fact can be found also for two-matrix models (see in this sense 
[p^]). But we dare say this problem is still open. 



5 The four— field KP hierarchy 

In this section we show how to compute the CF's of a definite model, the Ai'^, and 
of its reductions M\ and A^S], in genus by explicit use of the integrable hierarchy. 
As we shall see, once one knows the hierarchy of flow equations the CF's are almost 
completely known. The role of the H^-constraints is very limited. Let us start with 



3- 

5.1 The model Mj 

The model is specified by the Lax operator 

which arises in the context of the previous section when we consider a potential V2 of 
highest order 3. We have therefore altogether four fields 01,02,51,52. The Poisson 
brackets, the hamiltonians and some of the fiow equations have been given in [^]. 

The calculus of CF's via the hierarchy can be subdivided in several steps. 

Step 1. Meaning of the fields. We specify the connection of the fields character- 
izing the model with the correlation functions. Due to equation (2.16) we have 



ai =< rm > (5.2) 

In the following we will often denote ti by P. 

The interpretation for 02 can be inferred from the second flow equation for oi, 



^ = (a'i + 2a2 + 2ai5i)' (5.3) 



Therefore we have 



02 = ^ < T2P » < PPP » - < PP » 5i (5.4) 



via formal integration of eq.(5^). Formal integration is an operation we will often 
use in the following. It is the formal inverse of d (no integration constant involved) 
and is a proper operation in the context of integrable hierarchies whenever the pseu- 
dodifferential calculus applies. For example it is certainly correct when, as above, we 
apply it to expressions containing only the abstract symbols of the fields, it may be 
incorrect if applied to a particular solution. In the latter case one must be careful 
about the constants of integration. 
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Going back to eq. (|5.4| ), we see that the meaning of a2 in terms of correlation 
functions is clear once the meaning of Si is clear. An interpretation for the latter 
and for 5*2 can be derived from the first flows. In A4'^ the first flows are 

a'Q{N) = ai{N + l)-ai{N) 

a[{N) = a2iN + I) - a2iN) + ai{N)(ao{N) - ao{N - 1)^ 
a'^{N) = a2{N)[ao{N) - ao{N - 2) 

Let us introduce the operator Dq = e^" — 1 defined by 

DoJn = In+i - In 

for any function / depending on the discrete index N. Using the definition of the S 
fields, the first flows can be rewritten as 

S[ + DoS[ = Doai (5.6a) 
a'l = L>o«2 + aiDoSi (5.6b) 
a'2 = a2[DoSi + D0S2) (5.6c) 



It is very convenient at this point to consider to = N on the same footing as the 
couplings tfc and to associate to it an "operator" Q, in the same way as ti is associated 
to P, etc.. The label used above was motivated by the fact that do 



dto 



mimics 



the derivative with respect to A^. Keeping track of the operator Q is tantamount to 



studying the dependence in the CF's. Then (5.6a) implies 



5i =< [1 - e-^^j P >, ^2 =< e- 
Similarly, eqs.(5.6c) and ( 5.6b| ) become 



a2 



exp( «(e«-l)(l-e-"'^)» )/ 



2Q^ 



Doa2 =< PPP » - < PP »< (e^ - 1)(1 



^)P» 



(5.7) 

(5.8a) 
(5.8b) 



where / is an integration constant which in general might depend on all the couplings 
except ti. These two are compatibility conditions which come from the request that 
d and do commute. 

Eqs.(^.2|, ^^ ) and (5.7) provide the connection between fields and CF's we were 
looking for. 

Step 2. The genus contribution. In this section we will be interested in the genus 
part of the CF's, which corresponds to considering the dispersionless limit of the 
hierarchy, i.e. the limit in which, according to the degree analysis in section 3, only 
the first derivative terms (in ti and to) are retained (see |35],|Q). This corresponds 
to assigning the following degrees 



[tk] = x-k 
as in section 3, and moreover 



k>0, 



,{0)i 



.(o)i 



3, 



[SI 



(0)1 ^ [^(0) 



] = 1 
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and keeping the leading order terms in every equation. In this Hmit the two fields 
Si''^ and 82^^ collapse to a sing le field for 

5^ =« PQ »= 5f 

From now on in this subsection we drop the label which indicates the genus 
contribution since only the latter will be considered. With this understanding the first 
few dispersionless fiow equations are collected in Appendix Bl. The compatibility 
equations ( ^.8a ) and ( 5.8b| ) become 

a2 = e2«'3«»/o(to) (5.9a) 
< PPP »= doa2+ < PP »< PQQ » (5.9b) 



Moreover, from (|5.4| ), 

02 = ^ < T2P » - < PP PQ » (5.10) 

Step 3. The first critical point. For the type of solutions we are looking for the 
first critical point is by definition the one for which the correlation functions contain 
only nonnegative integral powers of the couplings. In this model this circumstance 
is implemented if 

3t3 = -1, tk = Q k>3 (5.11) 

We want to preserve homogeneity, therefore in the above degree assignment we set 
X = 3, i.e. 

[tk] = 3-k k>0, [ai] = 2, [02] = 3, [S] = 1, = 6 (5.12) 

The CF's will therefore be homogeneous functions of N, ti,t2- At times we will refer 
to this subset of the parameter space as the "small phase space" . 

This choice of the critical point corresponds, in the path integral language, to 
considering the potential 

Vl=tiTi+t2T2-^T3 (5.13) 

as the interacting potential, while the other are regarded as external fields. 
Step 4- Integrating the flow equations. From the homogeneity ansatz we have 

S = at2, ai = bti + ct2, a2 = dto + etit2 + ft2 

where a, b, c, d, e, / are numerical constants to be determined. If we plug these ex- 
pressions in the t2 dispersionless fiow equations 



^ = (2a2 + 2ai5)' (5.14) 

1^ = 245 + 4025' (5.15) 

Ot2 

f = (S^ + 2a,)' (5.16) 
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we find 



a = 2b, c = 2h\ e = 0, / = 

Then we can write 

< PP >= ]^{ati + a^tl) 

< T2P >= 2dto + a^tit2 + 

< PQ >= at2 

where, as in section 3, the symbol < ■ > denotes the correlation function evaluated 
in the small phase space. The second equation is obtained from ( ^.14 ) via formal 



integration. Integrating these equations with respect to ti,t2 and to respectively, and 
comparing the resultsQ, we find 2d = a and 

1 11 
< P>= -atj + atot2 + -ahitl + -a^t^ (5.17) 

After this example it should be clear how to proceed in order to calculate < T2>- 
We need to know < T2P >, < T2Q > and < T2T2 >. The first two are obtained by 
formally integrating the t2 fiows of oi and S. The third is obtained in the following 
way. We differentiate eq.( 5.1C| ) with respect to t2 and get 



da2 I ^ dS dai ^ 

dt2 2 dt2 dt2 



Then we plug eqs.( ^.14 ) and (|5.16| ) into this equation and obtain 



^ <r2T2P»= (4a25 + 2aiS2 + a2y (5.18) 

This can be formally integrated. So altogether we obtain 

< r2r2 >= 8025" + 4015"^ + 2al 

< T2P »= 2a2 + 2aiS 

< T2Q »= 52 + 2ai 

Evaluating the RHS's in the small phase space, integrating and comparing as before, 
we find 

1 1 

< T2 >= atoti + 20^*0*2 + a^titl + -a^tlt2 + -aHl (5.19) 

We can continue in this way and calculate the other one-point CF's, the pattern for 
the derivation does not change. They all depend on the constant a. This constant 
cannot be determined on the basis of the flow equations alone. Before passing to the 



*The reason why we have to compare the results of the different integrations, i.e. keep track of the 
integration constants, is that the correlation functions we find do not belong to the space in which is 
a formal integral (for example the space of rapidly decreasing functions). This remark applies to all the 
correlation functions calculated in this paper and means that integrability is effective in a far larger space 
than the one in which the formal rules of the pseudodifferential calculus strictly apply. 
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determination of o, let us remark that while the RHS of eq.( ^.15 ) is not a derivative 



w.r.t. ti, the RHS of (5.18) is. This circumstance, crucial for us to be able to calculate 



the CF (otherwise we would be left with an undetermined integration constant), is 
typical of all the flows and is rooted in the integrability of the hierarchy. 

Step 5. The W -constraints. To determine a we need the VF-constraints. Let us 
consider the W^lz^ = constraint at c = 0. It takes the form 

- <T2> +2t2 < n > +totl = 

Inserting the above expressions for < T2 > and < ti > we find that the constant a 
must be equal to 1. This is a general fact: provided a = 1, the correlation functions 
calculated by integrating the flow equations satisfy all the -constraints if c = 
(and the t2^k couplings are forgotten). This confirms what we have already found in 
section 3 concerning the NLS correlation functions. 

A sample of the one-point CF's is collected in Appendix C2. We recall that multi- 
point correlation functions involving only the operators Q, -P, T2 can be obtained from 
these results by simply differentiating a suitable number of times < P > and < T2 > 
w.r.t. to, ti and t2- For the multi-point CF's involving with k > 3, we have to first 
differentiate suitably chosen flow equations w.r.t. suitably chosen flow parameters, 
insert the flow equations in the expressions so obtained and proceed as above. 

A remark is in order concerning the constant a determined in this section. The 
value of this constant a = 1 is directly connected with the chosen fixed point Sts = — 1. 
The correlation function with generic a are appropriate to the choice 3^3 = — ^- In 
other words a is a normalization constant []. Nevertheless the fact remains that once 
the critical point is fixed the hierarchy is not enough to completely determine the 
CF's. We need also the 11^-constraints. 

Step 6. The CF < Q >. This is the only one-point function we have not 

calculated yet. From the above results we have 

<PQ>=t2, <T2Q>=ti+2tl 

Integrating the first equation w.r.t ti and the second w.r.t. t2 and comparing, we 
find 

<Q>= tit2 + ^ti + yto (5.20) 

where y is an undetermined numerical constant. This of course implies in particular 
that < QQ >= y. Further information could a priori come from the compatibility 
eqs.( ^.9a , 5.9b| ). However, inserting < PPP >= 1/2 and < PQQ >= 0, they only 



give us 

1 



/o = ^e~^'to 

Since we have used all the information at hand, the number y remains undetermined. 



t A normalization problem is always present in this kind of models and, when comparing CF's obtained 
in different contexts or with different methods, one should be careful about the normalization of both the 
coupling constants and the critical point. 
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Finally, integrating < Q >, < P > and < T2 > w.r.t. to,ti and t2 respectively, 
and comparing the results we find 



This completes our genus analysis of the model A4^. 



5.2 The model Ml- 

The model A4^ is obtained from A^l via hamiltonian reduction - the constraint is 
Si = -, and it is specified by the Lax operator 

L = + ai + a2Tr^ (5.21) 
a — D2 

We have therefore three fields ai,a2 and S2 = S. The Poisson brackets, the hamilto- 
nians and some of the flow equations have been given in [^]. 

The connection of the fields ai and S2 with the correlation functions is the same 
as before . The interpretation for 02 is slightly different from the previous one. In 
fact from the second flow equation for ai, |18|: 

dt2 ^ 

we deduce via formal integration 

02 = ^ < T2P > (5.22) 

The first critical point and the degree assignment are the same as in the previous 
model. 

Prom now on we will be dealing with genus and before we proceed let us sum- 
marize the basic relations. We have 



3*3 = -1, tk = k>3 

and 

[tk]=3-k k>0, [ai]=2, [02] =3, [S] = l, [J-(°)] = 6 
Moreover 

ai =< PP >, «2 = ^ < r2P », 5 PQ > (5.23) 

and again Q is coupled to = N. The CP's will therefore be homogeneous functions 
of N, ti, ^2. We have remarked above that this choice of the critical point corresponds 
to considering (|5.13| ) as the interacting potential. 

As before we now integrate the dispersionless flow equations (the first few are 
collected in Appendix B2). Prom the homogeneity ansatz we have as before 

S = at2, ai = bti + ct\, 02 = dt^ + etit2 + ft\ 
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where o, b, c, d, e, f are constants to be determined. If we plug these expressions in 
the t2 dispersionless flow equations 

dai , da2 / „\' dS / 2 , ^' 



we find 

a = b, ,c = e = 0, / = 

Then we can write 

< pp >= ati 

< T2P >= 2dto 

< PQ >= at2 

Integrating these equations with respect to ti,t2 and to respectively, and comparing 
the results, we find 2d = a and 

< P >= ^atj + atot2 (5.24) 

Proceeding now in the same way as before and, in particular, using the identifi- 
cation of a2 contained in ( 5.23| ), we find 



< T2 >= atoti + a^to^i 

as well as all the one-point functions we wish. All of them depend on the numerical 
constant a. 

It would be natural at this point to plug the one-point functions we found into 
the VF-constraints of subsection 2.2, as before. However one quickly realizes that 
this gives inconsistent results. Those W -constraints are incompatible with the re- 
duced hierarchy. This is not surprising. One should remember the origin of the 
M^-constraints. They come both from the string equations and the fiow equations. 
Since we have reduced the hierarchy it is evident that we have to change accordingly 
also the M^-constraints. Unfortunately we do not have a general theory for reduc- 
ing the W algebra and consequently the VF-constraints. However we can easily do 
without it. Remember that the only piece of information we need is the constant 
a. It is not difficult to guess the form of the constraints (at least in the case c = 0) 
which are compatible with the one-point CF's we have calculated. Once this form is 
determined it will tell us what the constant a is. In fact, if we consider the case c = 
and disregard the t2,k dependence, it is enough to calculate the Virasoro constraints, 

LnZ = 0|. 

The way to proceed is as follows. We try the following form for L_i 

1 00 o 

L-i = nT.^tkjz: + C-1 (5.25) 

2 dtk-2 



^The CF's of the model Ai^ do not satisfy Virasoro constraints such as LnZ — 0. Therefore some 
models satisfy Virasoro constraints, while others do not. This is a phenomenological observation we do 
not have an explanation for. 
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Here C_i is a term not containing derivatives w.r.t. the parameters. Now we impose 
L_iZ = 0, i.e. 



- <P> +2C_i = 

i.e. 

1 2 1 
C„i = jati + -atot2 

In a similar way we proceed with 

and with the other generators. We quickly realize that there exists a solution only if 

a = 2 



What we have done is legitimate if the new generators close over an algebra of 
the Virasoro type and if the constraints L^Z = are satisfied by the CF's for 
a = 2. These new constraints are then compatible with the reduced hierarchy. These 
requirements are actually true, to the extent we have been able to verify them. 
The new generators are written down in Appendix A2 and the first few correlation 
functions are collected in Appendix C3. 



5.3 The model M^. 

The model is obtained from via hamiltonian reduction - the constraint being 
S2 = -, and it is specified by the Lax operator 

L = d^ + aid + 02 (5.26) 

This is the Lax operator which gives rise to the Boussinesq hierarchy. We have 
therefore two fields ai,a2- The Poisson brackets, the hamiltonians and some of the 
flow equations have been given in In the Boussinesq hierarchy the flows 

with k = 1,2,3... do not appear |. It is therefore natural to ignore to = N. This 
is of course consistent only if we find a solution which does not depend on with 
k = mods. This remarkable fact indeed comes true. 

The connection of the fields oi and 02 with the correlation functions is the same 
as in the model 7W3 

ai =< PP >, «2 = ^ < r2P », (5.27) 

But the first critical point and the degree assignment are diff'erent. Precisely we find 
a solution with the desired properties if we set 

4*4 = -1, tk = > 4 

^From the reduction procedure from it is possible to define ^3,^6: •■■ flows in the reduced model, 
but they are not integrable. 
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and 



N=4-A: k>0, [ai]=2, H = 3, [^W] = 8 

We concentrate now on genus 0, as before. The CF's will therefore be homo- 
geneous functions of ti and t2- As before we now integrate the dispersionless flow 
equations (the first few can be derived from Appendix B3). From the homogeneity 
ansatz we have 

oi = at2, 02 = bti 

where a and b are numerical constants to be determined. Inserting this into the t2 
flow we flnd a = 2b. Using ( 5.27 ) we find 

< P >= atit2 



Similarly from 



we obtain 



Moreover 



2 

<T2P>=ati, < T2T2 >= --a^tl 



1 9 2 9 q 

< T2 >= -ati - -ahi 



1 9 9 2 q /I 

< r4 >= -ahit2 - —ahl 

and so on. As before we are left with an undetermined numerical constant a which 
only the appropriate PF-constraints can fix. 

We proceed as in the previous example to determine the VT-constraints appro- 
priate for this model. For example we guess the following form for c}-^^ 

-C-i = oE'^*fc7^ + C-i (5.28a) 

4" = oE'^^'^TiT+^o (5.28b) 

where the prime on the summation symbol means that the sum is over all > 
excluding the multiples of 3. Next we impose 

C^\Z = Q and 4^1z = 

We find that these conditions imply 

a = 6, C_i = 2iit2, Co = J 

A more thorough analysis confirms this and shows that the constraints compatible 
with the Boussinesq hierarchy are 

/:MZ = 0, r = l,2, n>-r (5.29) 
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\r] 

The generators Cn with r = 1,2 form a closed quadratic algebra, the W3 algebra, 
see Appendix A3. It corresponds to the version calculated in |3S] with central charge 
equal to 2. The first few one~point correlation functions of this model are collected 
in Appendix C4. 

Let us stress the distinctive features of the Boussinesq model: 

1) the CF's do not depend on A^; 

2) the appropriate W algebra closes over CH^ and d''; 

3) the coupling constant c does not come into play. 



6 Other models 
6.1 General recipe 

The three models discussed above describe essentially all the features of the A4p. 
And, as far as the method is concerned, it remains the same for all the other models. 
For all the A4p with 1 < / < p — 1 the first critical point is fixed by the condition 

ptp = -l, tk = k>p (6.1) 

and the degree assignment is 

[tk] =p-k, [J^]= 2p 

[S] = l, [ai] = 2,...,[ap_i] =p (6.2) 

where S denote the genus part of the fields Si, which is common to all of them. 
The CF's will be homogeneous functions oi to,ti, . . . , 

For the models A4p (the p-th KdV models) instead, one must first of all disregard 
all the tk with k a multiple of p; the first critical point is 

(p + l)tp+i = -1, tk = k>p+l (6.3) 

and the degree assignment is 

[tk] =p+l-k, [T] = 2p + 2 

[S] = l, [ai] = 2,...,[ap_i]=p (6.4) 

The CF's will be homogeneous functions of ti, . . . , tp_i. 

In all the cases the method consists in fixing the form of the fields on the basis 
of homogeneity; this leaves a few undetermined numerical constants; inserting the 
fields into the flow equations determines them up to an overall constant; in turn this 
can be fixed via the appropriate VK-constraints. The latter take the universal form 
(^.18 ) in all the unreduced M.^~^ models. In the other cases they can be constructed 



on the basis of the compatibility with the hierarchy characteristic of the model. 

6.2 Other examples 

For completeness we have explicitly worked out two more examples in genus 0, M\ 
and M.^. In the former the fields take the form 

2 12 2 1 

S = ts, ai = -t2 + 1% a2 = -ti + -^2*3 + 2*1, as = ^^o (6.5) 
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The first few dispersionless flow equations are collected in Appendix B4 and the first 
few one-point CF's are in Appendix C5. 
As for model A^4 we have 

ai = 12*3, 02 = 8t2, ag = 4ti + 18*1 (6-6) 

The first few dispersionless flows and one-point CF's are in Appendix B5 and C6, 
respectively. The generators for the appropriate VF-constraints are to be found in 
Appendix A4. 



7 Higher genus contributions 

Once the genus zero correlation functions of a given model are known, it is an easy 
task (at least in principle) to calculate the CF's in higher genus. There are several 
methods one can adopt. The most immediate one is to use the VF-constraints. These 
are either the original constraints ( p. 18 ) or the effective constraints we calculated in 



the previous sections. The important point is that, although the latter have been 
determined from the dispersionless versions of the relevant hierarchies, they are valid 
in general and provide the contributions to any genus. An example will suffice to 
illustrate this point. 

Let us consider the NLS model, which we met at the end of section 3 and was 
studied in detail in ||l^. Eq.( |3.10|) gives us the effective Virasoro constraints. We 
have seen that the genus zero contribution is represented by the highest order terms 



with respect to the degree assignment (|3.4| , p^ ). An inspection of eqs.( 3.10 , 3.11 ) 
tells us that the second nontrivial contribution will have degree 2x less than the 
leading order, the third 4x less, etc. Therefore the free energy has a genus expansion 

oo 

^=^^W, [^C^)] = 2x(l - /i) (7.1) 

h=0 



Inserting this into ( 3.10| ) we find 



n-l 



^ ktk < Tk+n »h +2N < X„ >/, + ^ ( < TkTn-k 
k=2 k=l 

h 

+ ^ < Tfc Tn-k :>h-h' ) + {N^Snfi + Nti6n,-l)Shfi = (7.2) 

h'=0 

where <^ • represents the genus h contribution to the <^ ■ » CF. Eqs.(7.2) 
are recursive, therefore it is easy to turn the crank and calculate the higher genus 
contributions 

to the order we wish from those of genus 0. A few examples at the critical point 
2*2 = — 1 are given in Appendix D. 

Prom the structure of the H^-constraints appearing in two-matrix models it is easy 
to see that genus recursiveness is a general characteristic. Therefore computability 
of higher genus contributions is guaranteed. We recall that genus recursiveness is 
typical of topological field theories. 

Clearly it is possible to calculate higher genus correlation functions starting from 
the relevant hierarchies and computing the genus expansion of them. This is the best 
method if one wishes compact formulas for the CF's. 
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8 Higher critical points 



The first critical point by definition implies a dependence of the basic fields on the cou- 
plings specified by homogeneous polynomials (i.e. with non-negative integer powers). 
Higher critical points are characterized still by a homogenous dependence, but with 
rational and/or negative powers of the couplings. For the models of type Mp with 
I = 1, . . . ,p — 1, characterized by a small phase space with parameters to,ti, tp_i, 
higher critical points are specified by 

ktk = ±1, k> p, ti = ioT l> k (8.1) 

and the degree assignment is [| 

[ti] =k-l, ,[^]='^k 

[S] = l, [ai] = 2,...,[ap_i] =p (8.2) 

For the models of type Aip (p-th KdV hierarchy) characterized by a small phase 
space with parameters ti, ...,tp_i, higher critical points are specified by 

ktk = ±1, k>p+l, ti = Q ioi l>k, k^np (8.3) 

and the degree assignment is f\ 

[ti] = k-l, ,[J^]=2k 

[ai] =2,..., [ap_i] = p (8.4) 



The lb sign in eqs.(8.1, |8.3| ) has to be chosen in such a way as to avoid complex 
coefficients in the small phase space expressions of the basic fields. 

A few examples will suffice. For the model A4^ the second critical point is fixed 

by 

5t5 = 1, ti = l> 5 
and, for simplicity, we set ^4 = 0. Consequently the degree assignment is 
[ti]=5-k, m = W, [ai]=2, [02] =3 

It is convenient at this point to take immediately into account the constraint C^^^Z = 
0. It reads 

< T2 >= -6*1*2 (8.5) 



Now eq.( 5.27 ) implies 02 = —3*2- Using again the second flow equations (Appendix 



B3) we also find 



1 



ai = 3tl (8.6) 



*The parameters to-,ti, . . . , ip_i are the self-coupling parameters of the model, while tp, . . . , tk-i are 
to be regarded as external couplings. If one wants models in which all ti, t2, ■ ■ ■ , tk~i are self-interacting 
parameters one should look at models of the type TWj.. 

'''The tp+i, . . . , tk-i are to be regarded as external couplings, see previous footnote. 
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At this point it is easy to calculate the correlation functions following the method of 
section 5. A few of them are collected in Appendix C4. 

Another example is the second critical point of the model M4: 

6*6 = -1, ^5 = 0, ti = l>6 
Consequently the degree assignment is 

[ti]=6-k, [;P] = 12, M = 2, [a2]=3, [03] = 4 
The constraint C^^\z = tells us that 

< T2 >= 8tl + 12*1*3 (8.7) 
Prom this and the second flow equations we can derive 

ai = ~ht^\ 02 = 6*3, 03 = ^*?*3"2 + 4*2 (8.8) 

from which we can easily derive the CP's (see Appendix C6). 

We conclude that we can at least in principle (i.e. modulo technical difficulties) 
solve the Mp models at higher critical points as well. 



8.1 Non— perturbative equations 

This paper privileges genus expanded solutions. Wc do not want to get involved here 
into the problems of non-perturbative solutions. However it is useful to spend a few 
words at least on the way we can extract non-perturbative equations. We show a 
few examples. The first is the NLS model {M2 model) at the second critical point 
3*3 = —1, *2 = 0. Differentiating the wj^lz = at c = w.r.t. *o and *i and using 
the second flow equations one finds 

R' + 2RS = to, -S' + 2R + S^ = ti (8.9) 

where, as usual, we set ai = R. These two equations must be supplemented with the 
*o flows: 

(e^°-l)i? = e^"5', (e^o - 1)5 = (Ini?)' (8.10) 

The above equations specify the *o,*i dependence of R and S. 

The second example is the second critical point of the M^, model, which has been 
specified above. Differentiating w.r.t. *i and *2 the constraint C^}^^ and using the 
second flow equations we obtain 

7 2 

2a2 - a[ + 6*2 = 0, 44 - -a'l - -aj + 6*1 = (8.11) 

o o 

These should be supplemented with the (complete) second flow equations 

^ = 24 - a'l, ^ = 4' - -a'C - -aia[ (8.12) 
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The four equations above specify the non-perturbative dependence of 01,02 on ti 
and t2. 

An interesting example is the next critical point in the same model, Al^j, i.e. 



7t7 = -1, ^4 = 0, ti = 0, l>7 (8.13) 

il CC 

to ti and t2 we obtain 



and is is an external coupling. Differentiating the constraint C^^iZ = with respect 



- < T4P > +5*5 < T2P » +6*2 = 

- < T4r2 > +5*5 < r2T2 > +6*1 = 

respectively. Now, using the *2 and *4 flows (see Appendix B3), we find 

a'l' - 2o2 + 2aia'i - 4aia2 + 15*5(2a2 - a[) + 18*2 = (8.14) 

1 8 12 

-a'l" - 402 + 2aiai + a[a[ + 4o'i02 + -af + 15*5( — o" - -aj) + 18*i = 

3 3 2 3 

With the redefinition 

15 1 , 

oi ^ oi - Y*5, 02 ^ 02 - -ai 

and a shift of *i, we obtain the same equations derived in ref.||3| (see also [^]) for the 
Ising model on discretized Riemann surfaces (up to normalization of the quantities 
involved): *i,*2 and *5 are related to the cosmological constant, the magnetic field 
and the parameter T of respectively. Eqs.( |8.14 ) together with ( 8.12 ) determine 
the *i,*2 dependence of 01,02. 



9 Interpretation and comments 

Let us briefiy introduce in this last section two points that require a longer and more 
careful elaboration: a possible connection of the previous models with 2D gravity 
and topological field theories. We will then comment on some leftover problems. 

On the basis of the continuum theory one would expect the chemical potential 
(to be assimilated to oi) of 2D gravity coupled to conformal matter in genus 0, to 
behave like 

2 

ai ~ *po+9o-i (9.1) 

where * is proportional to the cosmological constant, and po,qo are relatively prime 
integers. From the above examples and general statements one can see that the 
behaviour of the models Mp at the various critical points can abundantly account 
for (|9.1j) , provided we interpret either *i or *o as *. 

Let us come now to the interpretation of the models we have presented in this 
paper in terms of topological field theories coupled to topological gravity. The latter 
models are completely specified once one gives the primaries Oa, the metric ??q,^ and 
the fusion coefficients Caff . As primaries of our models we can take the fields coupled 
to the parameters of the small phase space, rj and C can be derived from the CF's of 
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the primaries at the first critical point. For example, for the model Ai]^ the primaries 
can be chosen to be: tq = Q, ti = P and T2 and by inspection of Appendix C3 we 
find 



V02 = Vn = mo = ^ 

while the remaining Tjap vanish. For some models, like the Aip ones, a similar iden- 
tification can be very easily carried out (see also |2^). There is therefore room for 
an interpretation in terms of topological field theories. However further elaboration 
is required. We defer a complete discussion on this point as well as a study of the 
connection with Landau-Ginsburg theories to a forthcoming paper. 

It remains for us to point out a few problems concerning the completeness of 
the proofs and arguments presented in this paper. The first remark concerns the 
statements, made in section 5 and 6, about the agreement between the calculated 
correlation functions and the Vl^-constraints. The agreement has been checked up to 
a remarkable order (seventh, eighth or even higher order CF) with a computer. In 
order to produce a rigorous proof one should be able to find compact formulae for 



the CF's, analogous to those we exhibited in |16]. These formulae are perhaps not 
beyond reach. 

A second question is: did we study in this paper all the models contained in 
the two-matrix theory? We recall that for any small phase space of the type 
ti,t2, ■ ■ ■ , tp-i or to, ti, . . . , tp-i, for any p, we presented one or more models. Does 
two-matrix theory describe other integrable models with a more peculiar small phase 
space (such as, for example, ti, ^2, ^4, tj)7 To answer this question one has to ascertain 



whether there exist other reductions beside those considered in |18],|21]. 

An entirely new problem is posed by the calculation of the correlation functions 
with c ^ 0, as we did in section 3, but via the integrable hierarchies. In this case 
we have to consider two integrable hierarchies, the first (system I) and the second 
(system II). The novelty consists in the fact that we want to know the dependence 
of the solutions of the first on the parameters of the second and vice versa, plus the 
dependence on c. In we showed that the corresponding fiows commute. 

Finally a few words on general multi-matrix models (with bilinear couplings). 
The results of [0] compared with those of |21] tell us that very little has to be 



expected in terms of new integrable hierarchies and new models. Of course if one 
wishes to study CF's dependent on all the couplings (see the previous paragraph) 
the problem presents new and perhaps interesting aspects. 

Appendices 
Appendix Al 



The basic >Cn (1) generators that appear in ( 2.1^ ) are 



£[■1(1) . j:«,,_^+(iv+!i^)-^+i2 



k=l 



+Nti,i5n,-i + ^N{N + l)5n,o, n>-l. 
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and 



4^1(1) ^ y: — + E^*M E 



h,/2=l 



3 ^ ^ dti^idti^kdti^n-i-k 

n—l o2 °° r) 

+(^ + ? + 1) E ^71^ + (2Ar + n + 2) ^ lU 



n 

2 dtkdtn-k dtn+i 



k=l 

,2 - . ..... . 1. ... ... d 



+iN' + {n + 2)N + -{n + l)(n + 2))- 

+(2NHr^2 + Ntl{)5n,-2 + N{N + l)ti,i<5„,_i 

+iiV(iV + l)(iV + 2Xo, n>-2. 

In the above formulas it is understood that when ti^k appear with /c < 0, the corre- 
sponding term is absent. As already noted the higher rank generators can be derived 
from the algebra W itself. 

Actually it remains for us to specify what the £{[''(l)'s are. These are given by 

4°] = ^ + iV5n,o, n > 

These generators represent an abelian extension of the W algebra, of which they form 
an abelian subalgebra. The other generators behave tensorially with respect to £n'. 
However since the latter play a minor role in this paper, we do not insist on this 
point. 

Appendix A2 

Here are the generators of the Virasoro algebra for the reduced model 

1 °° d 1 

= ^E^^.^n^? + ^*- 

" 2^^ Otfc 8 16' 



dtk+2n 4 dt2n 8 dtkdt2n-k 



Appendix A3 

The generators of the W algebra appropriate for the (Boussinesq) model are 

k+l=3n k+l=-3n 
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d 



^ <JHi+l2+2,n ^' 



92 



h,l2 = l 



"^243 ^ 



l—k—j=—3n 



dti^dtj 



dtidtkdtj 9 

l.k,j -' l,f^-j 

l+k+j=3n l+k+j=—3n 



Vn 



In these expressions summations are limited to the terms such that no index involved 
is either negative or multiple of 3. Contrary to the previous cases we have introduced 
the generators for all n. In this way the above two sets of generators form a closed 
algebra, the W3 algebra, 0,11], @ 





[2] 



n+m 



(2n — m)C 

— — (n — m)[{n^ + w? + Anm) + 3(n + m) + 2j£j^^ 
y oiu 



[1] 



where 



k<-l k>0 



This corresponds to the quantum W3 algebra calculated in [33| when the central 
charge is 2. 



Appendix A4 

The generators of the W algebra appropriate for the model are 

= \flktk^^ + ^T + ^ ^ kltkti + ^6n,o, Vn 

4 fc^i dtk+4n 32 dtkdUn-k 2 32 



l+k+j=8 

In these expressions summations are limited to the terms such that no index involved 
is either negative or multiple of 4. 

Prom the above generators we can generate the W4 algebra. 
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Appendix Bl 

Here are the first few equations of motion of the model Ai^ in the dispersionless hmit 
(tk = h,k)' 



dai 

ds_ 

dt3 



(SaiS^ + 6a2S + Safj' 
3a2>S'^ + 602(5'^)' + 3aia2 + 60102 
(5^ + 6015 + 302)' 



dai 
da2 

ds_ 



(4oi5^ + 120252 + 12o?5 + 120102)' 
: Aa'^S^ + 802(5^)' + 12010^5 + 24o2(oi5)' 
(S* + 12ai5^ + 12025 + 6a?)' 



Appendix B2 

Here are the first few equations of motion of the model M\ in the dispersionless limit 
and with the convention = ii,fc: 

— = 2o'2, ^=2(a25), ^ = {s^ + ai) 



|^ = (^a? + 3o25)' 

|g=(^Oi02 + 30252)' 



dai 
'dtl 

902 

dS_ 

5*4 



: (4025^ + 40102)' 
: (402 5^ + 2o^ + 4oi025)' 
(5^ + 2oi52 + 4o25 + o?)' 
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Appendix B3 

Here are the first few (complete) equations of motion of the model with the 
convention tk = ti,jfc. We disregard the tk flows with k = mod 3: 

dai I ,1 da2 „ 2 ,„ 2 , 
— = 2a,-a„ gf^=a,-^a, - -a,a, 

dai / I ,„ 2 „ 2 , A \i 

da2 /I /// 2 ,,,, 2 / 2 2 2 „ 1 / / ^ 3V 

Appendix B4 

Here arc the first few equations of motion of the model AI4 in the dispersionless limit 
{tk = tl,k)'- 

^ = (2a2 + 2aiS)', ^ = 2a^ + 2^5 + Aa2S', 

dt2 ^ ^ dt2 



(Saj + 3a3 + SaiS'^ + 602 5)' 



3a2S^ + 602 (5^)' + Ga'^S + Dag-S' + Sai^ + 6a;a2 



3aia3 + 9aia3 + BagS^ + 903(5^)' 



H = (53 + 6015 + 802)' (9.2) 

Appendix B5 

Here are the first few equations of motion of the model M\ in the dispersionless limit 

and with the convention tk = ti,fe- We disregard the flows with = mod 4: 

dai „ / da2 /„ 1 2\' ^^3 1 , 

W2 = m-2 = y^""' - 2^0 ' or2 = -2"^^^ 

dai / 3 2,0 y 3 , 5a3 3 2n/ 3 , 

— = (--«, + 3a3) , ^ = --(aia2) , — = --(a2) + -a^as 



— = (-a2 + -aia, - -a,) 



da2 /5 15 2 

5 2x/ , 5 2 / 5 , 
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Appendix CI 

In this Appendix we collect the first few genus one-point CF's calculated in section 
3. Small phase space: N,ti. Critical point: 2t2 = —1. 



< Tl > 

<T2> 
<T3> 

< r4 > 



N 



t + cs 



<ai>=N 



s + ct 



s2 + t2 + 2cst N"^ - Ns^ 
N — + 



N- 



(1 - c2)2 1 - c2 

^^^^ ((?s^ + 2cst + t^ + m - SATc^) 



(1 - C2) 

N 



c^s'^ + Ach^t + Qc^s^t^ + Acst^ + 



(1 - C2)4 

+mc^s^ - QNc*s^ + UNcst - UNc^st 
+6Nt^ - mcH'^ + 27V2 - 4iV^c2 + 2N'^c^') 



and in general 



Appendix C2 



< CTfc > (i,s) =< Tfc > (s,t) 



One-point correlation functions of the model in genus 0. Small phase space: 
to = N,ti,t2- Critical point: 3*3 = —1. 



<Q> = tit2 + -4 + yto, 



<p> 



-tl + tot2 + 2*1*2 + 



<r3> = \4 + ^toht2 + ^tl4 + ^h4 + ^4 + Ato4 + 4, 



<n> = ^totltlt2 + 9toti4 + ^ti4 + ^ti4 + 3*0*2 + y + 24> 



^ 2 2 3 3 2 2 4 

< Ts > = -igti + —tQitit2 + 25*0*1*2 + ^*1*2 + y *1*2 

5/1 45 fi , „ 9 9 37 c 65 o 

+ ^t\ + -*i*6 + 10*2*2 + —to4 + -^8 



In the first equation y represents an undetermined numerical constant. 
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Appendix C3 

One-point correlation functions of the model M.^ in genus 0. Small phase space: 
to = N,ti,t2- First critical point: 8*3 = —1. 



<Q> = 2tit2 + -4 + yto, <P> = tl + 2tot2, <T2> = 2toti + 4:to4 



<T3> = 4 + 6totlt2 + 8to4 + 7^0 



<T4> = 4tot? + 16*0*1*1 + 4*g*2 + 16*0*2 

< 7-5 > = ^*i + ^*i*o + 15*o*?*2 + 40*o*i*i + 15*g*i + 32*0*^ 

Appendix C4 

One-point correlation functions of the Boussinesq model in genus 0. Small phase 
space: *i,*2. First critical point: 4*4 = —1. 

< n > = 6*1*2, <T2> = 3*? - 8*i 

(9.3) 

<T4> = 12*^*2 - 16*^, <T5> = 5*f - 40*i*i 

256 

<T7> = 28*?*2 - 112*1*^, <T8> = lOtf - 160tl4 + 

Here follow the one-point correlation functions of the same model at the second 
critical point: 5*5 = 1. 



3 

< n > = 2tl - 3*1, <T2> = - 6*1*2 



3 5 
<T4> = -8*1^*2 + 4*1, <T5> = - 2tl + 15*1*1 



(9.4) 



Appendix C5 

Genus one-point correlation functions of the model M.\. Small phase space: *o 
A?^, *i,*2,*3. First critical point: 4*4 = —1. 

<Q> = *i*3 + ^ii + 3*2*1 + 7*1 + y*o 



2 , 4ol0o„. 

<P> = -*1*2 + *l*i + *0*3 + 3*1*3 + ^^2*1 + 2*1 
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1 2 8 10 , 4 o 2 8 , 

<T2> = -tt + ^ht2t3 + Yt^4 + -tot2 + 3totl + —tl 

16 2 2 34 4 19 e 



<T3> = 10*1*2*1 + 10*1*3 + *1*3 + toh + ilti + 6*0*2*3 

+9*0*1 + y *ii3 + f *2*i + 38*2*1 + 19*^ 

Appendix C6 

Genus one-point correlation functions of the model M^. Small phase space: 
*i,*2,*3- First critical point: 6*5 = —1. 

< n > = 8*^ + 12*1*3, <T2> = 16*1*2 - 72*2*1 

<T3> = 6tl- 72*^*3 + 81*1 



< rs > = 40*1*1 + 30*^*3 - 360*^*1 + 243*| 

One-point correlation functions of the same model M4 at the second critical point: 
6*6 = -1. 

< n > = 12*2*3 - hlt^\ <T2> = 12*1*3 + 8*2 

o 

<T3> = 12*1*2 + ^*?*3 2 - y*i 
<r5> = -^*t*3'-y*?*2*3'' + Y*l*l + 30*1*3 

9.1 Appendix D 

In this Appendix we collect results on correlation functions of the NLS {M2) model. 
Small phase space: iV, *i. First critical point: 2*2 = —1. The genus h contribution is 
denoted < ■ >h- 
Genus 

2<2ifc<r+2 

Genus 1 



< Tl >1 = < T2 >1 = < T3 >= 0, < TA >= N, 
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< r5 >i= 5iVii, < re >i= IbNtj + lOiV^ 

< T7 >i= 35Ntl + 70N\, < T8 >i= lONtf + 280N^tl + 70N^ 

< rg >i= 126Ntl + 840iV2tf + GSONHi, 

< no >i= 210Artf + 2100Ar2tf + 3150iV^t? + 420Ar^ 

(9.5) 

< nn >i= 0, 1 < z < 4, < nrs >i= 5N, < nre >i= 30Nti 

< T2T2 >i=< r2r3 >i= 0, < T2T4 >i= 4N 

< T2T5 >i= SOATti, < T3T3 >i= 3iV 

Genus 2 

< n >2= l<i<7, <Ts >2= 21iV, < rg >2= 189iVii 
< rio >2= 945iVt? + 483Ar2, < m >2= 3465Art? + 5313Ar2ti 
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